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■ Aims. We study the hard X-ray luminosity function and absorption distribution of local (z < 0.1) active galactic nuclei (AGN) 
and discuss the implications for AGN cosmological evolution and for the cosmic X-ray background (CXB). 
£X| Methods. We use the INTEGRAL all-sky hard X-ray survey to perform a statistical study of a representative sample of nearby 
AGN. Our entire all-sky sample consists of 127 AGN, of which 91 are confidently detected (> 5a) on the time-averaged map 
OO obtained with the IBIS/ISGRI instrument and 36 are detected only during single observations. Among the former there are 66 
non-blazar AGN located at \b\ > 5°, where the survey's identification completeness is ~ 93%, which we use for calculating the 
' AGN luminosity function and X-ray absorption distribution. 
Results. In broad agreement with previous studies, we find that the fraction / a of obscured (log Ah > 22) objects is much higher 
(~ 70%) among the low-luminosity AGN (Lhx < 10 43 ' 6 erg s _1 ) than among the high-luminosity ones (Lhx > 10 43 ' 6 erg s _1 ), 
/ a ~ 25%, where Lhx is the luminosity in the 17-60 keV energy band. We also find that locally the fraction of Compton-thick 
AGN is less than 20% unless there is a significant population of AGN that are so strongly obscured that their observed hard 
i-S^ • X-ray luminosities fall below ~ 10 40 -10 41 erg s _1 , the effective limit of our survey. The constructed hard X-ray luminosity 
Or function has a canonical, smoothly broken power-law shape in the range 40 < log Lh x < 45.5 with a characteristic luminosity of 
q log L* = 43.40 ± 0.28. The estimated local luminosity density due to AGN with log L hx > 40 is (1.4 ± 0.3) x 10 39 erg s" 1 Mpc~ 3 
!— i ' (17-60 keV). We demonstrate that the spectral shape and amplitude of the CXB are consistent with the simple scenario in 
which the Nn distribution of AGN (for a given L^/L t (z) ratio) has not changed significantly since z ~ 1.5, while the AGN 
Cd luminosity function has experienced pure luminosity evolution. 
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1. Introduction obscured sources. Furthemore, since deep extragalatic 

X-ray surveys usually cover small areas on the sky, they 
Recent deep. extragalatic X-ray surveys (see do nQt pr£)be the bfight end of ^ AQN luminosity 

i Brandt fc Hasinger | |2005| for a review) have signifi- distribution at low redshift {z ^ 3) It is m these 

cantly advanced our understanding of the cosmological twQ overlapping areag _ heavily obscured AGN and 

evolution of active galactic nuclei (AGN) and the growth luminous nearby AGN _ where deep x _ ray surveys could 

of massive black holes. In particular, these surveys have be efficiently comp i em ented by all-sky hard X-ray (above 

revealed the phenomenon of AGN "downsizing" , i.e. ^ iq keV) surveys 

gradual transition from the dominance of luminous 

quasars at high redshift (z £ 2) to that of lower lumi- Recently, a serendipitous all-sky survey in the 3- 

nosity AGN at z & 1. Another important result is the 20 keV energy band was performed with the RXTE obser- 

discovery of numerous type 2 (i.e. with significant X-ray vator y ( the RXTE Slew Survey, or XSS), in which ~ 100 

absorption) AGN at medium and high redshifts. Despite AGN were detected at |6| > 10° (|R.evnivtsev et al.||2004|). 

these successes, the census of AGN in the Universe may This sample was used to construct an X-ray luminosity 

still be significantly incomplete, in particular because function of local (z < 0.1) AGN and to study their absorp- 

deep X-ray surveys, performed effectively at energies tion distribution (|Sazonov k Revnivtsev||2rjOJ). However, 

below - 8 keV, are biased against detecting heavily despite its harder X-ray range compared to traditional 

X-ray surveys, the XSS was still biased against finding 

Send offprint requests to: sazonov@mpa-garching.mpg.de sources with absorption columns Nh > 10 23 cm~ 2 . 
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A comparable or better (for heavily obscured sources 
and/or in crowded sky regions) capability for studying 
bright (j> 1 mCrab) hard X-ray sources is no w pro- 
vided by the IBIS telesco pe JlJbertini etal]l2003l) of the 
INTEGRAL observatory l|Winkler et alJl2003|) . owing to 
its good sensitivity above 20 keV, large field of view 
(~ 28° x 28°), and good (~ 10') angular resolution. Since 
its launch in October 2002, INTEGRAL had observed 
~ 60% of the sky by the spring of 2005, mostly around 
the Galactic plane and the Galactic center. This moti- 
vated us to initiate a series of pointed INTEGRAL obser- 
vations (Proposal 0320108) aimed at covering the remain- 
ing "empty" (extragalactic) sky regions. This campaign 
has now been completed and practically full coverage of 
the sky has been achieved after all public and our propri- 
etary INTEGRAL observations were added. 

Based on the all-sky hard X-ray map obtained and the 
follow-up source identifica tions, we co mpiled a catalog of 
detected sources IjKrivonos et al.ll2006l) . In the present pa- 
per we use the INTEGRAL all-sky survey to study such 
key statistical properties of the local AGN population as 
the hard X-ray luminosity function and absorption distri- 
bution. We also discuss the implications for the cosmo- 
logical evolution of AGN and for the cosmic X-ray back- 
ground. We note that there have been previous attempts 
to use INTEGRAL fo r studying the statistics and demog- 
raphy of nearby AGN jKrivonos et alJ20o"5llBassani et alJ 
l200alBeckmann et al.l2006|) . The significant advantage of 
our present study is that it is based on an all-sky, largely 
serendipitious survey, which allowed us to significantly 
alleviate the previously significant bias associated with 
INTEGRAL AGN pointings. We also note that a similar 
all-sky hard X-ray survey has since recently been car ried 
out by the Swift observatory l|Markwardt et alJl2005h . 

2. The AGN sample 

The catalog of iKrivonos et all l)200(j[) includes in total 127 
AGN, which are listed in Tabled Of these, 91 AGN (re- 
ferred to below as the main sample) were detected with 
more than 5ct significance on the time-averaged all-sky 
map obtained with IBIS/ISGRI in the 17-60 keV en- 
ergy band and are thus well-suited for statistical analy- 
sis. In addition, there are 36 AGN, listed in the lower 
part of Table ^ that were confidently detected in single 
INTEGRAL orbits or observations (~l-3 days long) but 
fall below the 5cx threshold on the average map. We ex- 
clude these transiently detected AGN from the following 
statistical analysis. 

Most of the AGN presented were known as such before 
INTEGRAL, while 33 sources (those with "IGR" names) 
have been identified as AGN primarily based on their hard 
X-ray detection with INTEGRAL. Two of these newly dis- 
covered AGN - IGR J09446-2636 and IGR J13149+4422 
- are reported here for the first time. As indicated in col- 
umn 2 of Tabled] 77 objects in our main sample are opti- 
cally classified as Seyfert galaxies (4 of them radio loud) 
and 9 are blazars (BL Lac objects or flat-spectrum ra- 
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Fig. 1. Hard X-ray luminosity vs. redshift for the identi- 
fied AGN from the INTEGRAL all-sky survey (the main 
sample). Circles and squares indicate emission- line AGN 
and blazars, respectively. Empty symbols indicate AGN 
located at |6| < 5°, which are not used in the statistical 
analysis. 

dio quasars). Although the optical class has not yet been 
reported for 5 objects, their X-ray, optical, infrared, and 
radio properties strongly suggest that they are all Seyfert 
galaxies. The quoted optical classes are based on the NED 
unless a reference is provided. Of the AGN at |6| > 10° in 
our main sample, 30 were also detected during the RXTE 
Slew Survey. 

We found published rcdshifts for all but 4 of our AGN 
(column 3 of Table and used these to determine the 
source luminosities in the observed 17-60 keV energy 
band (column 7) from the hard X-ray fluxes measured 
by IBIS/ISGRI (column 6). For the nearest sources (z < 
0.01) the di stances from the Nearby Galaxies Catalogue 
l|Tiillvlll988h were adopted (column 4) , while the luminos- 
ity distances to the other objects were determined from 
their redshifts assuming a cosmology with fl m = 0.3, 
= 0.7, and H a = 75 km s _1 Mpc -1 , which is 
used throughout the paper. Figure ^ shows the distribu- 
tion of our AGN (of the main sample) in the redshift- 
luminosity plane. One can see that all of the emission-line 
(i.e. non-blazar) AGN are located at z < 0.1, except for 
IGR J09446-2636 at z w 0.14, while there are several 
blazars as distant as z ~l-2.5, with (isotropical) lumi- 
nosities up to ~ 10 48 erg s _1 . Because of the proximity of 
the non-blazar AGN, our quoted values for their observed 
luminosities may also be regarded as their rest-frame 17- 
60 keV luminosities. 
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Table 1. INTEGRAL all-sky survey AGN sample 



3 



Object Class" z D Note 5 -F17-60 kov log L hx Nh Ref° 



(sources at \b\ < 5° are given in bold) Mpc 10 11 crg/s/cm 2 crg/s 10 22 cm 

AGN detected at > 5cr on the time-averaged IBIS/ISGRI map 



Mrk 348 


S2 


0.0150 






7.41 ±0.83 


43.52 


30 


2 


IGR J01528-0326=MCG -01-05-047 


S2 


0. 


.0172 




1 


1.63 ± 


.30 


42.98 






NGC 788 


S2 


0. 


0136 






4.83 ± 


.29 


43.25 


44 


3 


LEDA 138501 


SI 


0. 


.0492 






3.96 ± 


.72 


44.30 


< 1 


4 


Mrk 1040 


SI. 5 


0. 


0167 






4.85 ± 0. 


.81 


43.43 


< 1 


5 


IGR J02343+3229=NGC 973 


S2 


0. 


.0162 




1 


3.88 ± 0. 


.63 


43.31 






NGC 1068 


S2 


0. 


0038 


14.4 




1.90 ± 0. 


30 


41.68 


> 10 3 


6 


ATT it *~i A 1 | £} -t 

4U U241+bl 


SI 


0. 


.0440 






4.75 ± 0. 


.60 


44.28 


< 1 


7 


IN LiO 1142 


S2 


0. 


0288 






4.61 ± 0. 


.35 


43.89 


45 


8 


Iri Uozo+o4z 


SI 


0. 


0610 






2.74 ± 


.47 


44.34 


< 1 


4 


TS.T /■""' /"i 1 o c r 
NGC 1365 


Si. 8 


0. 


0055 


16.9 




3.30 ± 


.65 


42.06 


~ 50 


9 


3C 111 


BLRG 


0. 


.0485 






7.83 ± 0. 


.87 


44.59 


< 1 


10 


TT 1 /~i noo /~< nno 
LibO 033-G002 


S2 


0. 


0181 






1.94 ± 


.27 


43. 10 


1 


1 1 


IKAb 05078+ lb2b 


Si. 5 


0. 


0179 






5.93 ± 


.75 


43.58 


< 1 


4 


Mrk 3 


S2 


0. 


0135 






6.83 ± 0. 


.30 


43.39 


110 


6 


Mrk 6 


SI. 5 


0. 


0188 






3.66 ± 


.30 


43.41 


~ 5 


12 


IbU J07503— 4137= 2MAbA J075bl9b3— 4137420 


S2 


0. 


0210 




2,3 


1.24 ± 


.24 


43.04 


< 1 


13 


IGR JU7597— 3842 


Si. 2 


0. 


0400 




3 


2.91 ± 


.26 


43.98 


< 1 


8 


LbO 209-0.012 


SI. 5 


0. 


0405 






1.66 ± 


.24 


43.75 


< 1 


4 


TT 1 ^ * „ 11 A A A G 

rairall 114b 


Si. 5 


0. 


0316 






1.62 ± 


.25 


43.52 


< 1 


4 


IKAb 09149 — 6206 


SI 


0. 


0573 






2.06 ± 


.26 


44. 16 


< 1 


4 


A ^" ™1 1 1 n 

Mrk 110 


NLS1 


0. 


0353 






5.86 ± 1. 


. 14 


44. 17 


< 1 


3 


Ti^i o TnnH ncoc c J "rnn^i/ioTn o/^oor/^ 

IGR J 0944b — 2bob— bar 1 J 0944370 — 263356 


SI 


0. 


.1425 




4 


3.91 ± 


.74 


45.27 


< 1 


4 


NGC 2992 


S2 


0. 


0077 


30.5 




5.08 ± 0. 


.39 


42.76 


1 


14 


A .T z"' c no i n 

MOb -5-23-lb 


S2 


0. 


0085 






9.77 ± 0. 


.84 


43. 14 


2.3 


15 


TS.T /■""< /"i o n o i 

NGC 3081 


S2 


0. 


0080 


32.5 




4.62 ± 


.56 


42.77 


50 


14 


boO 203-G013 


S2 


0. 


0333 






2.15 ± 


.39 


43.69 


40 


16 


TS.T THOT 
NGC 3227 


Si. 5 


0. 


0039 


20.6 




9.05 ± 0. 


.84 


42.67 


< 1 


3 


tvt 1^1 o n O 1 

NGC 3251 


S2 


0. 


0107 






3.87 ± 


.64 


42.94 


150 


17 


Ti^i t> t i no oc Art A T om A ov iiaoojcoa ^rurcoi 

IGR J 1038b — 4947— 2MAda J 10384520— 494b531 


SI 


0. 


.0600 




5 


1.47 ± 0. 


.24 


44.05 


1 


8 


IGR J10404— 4b25— LEDA 93974 


S2 


0. 


0239 




6 


2.11 ± 0. 


.35 


43.38 


2.8 


8 


TVT (~S (~\ OI7DQ 

IN GO 3783 


SI 


0. 


0097 


38.5 




12.28 ± 1 


.86 


43.34 


< 1 


3 


IGR J 1202b — 5349— WKK 05b0 


S2 


0. 


.0280 




2,3 


2.44 ± 


30 


43.59 


2 


13 


NGC 4151 


SI. 5 


0. 


0033 


20.3 




■IT o o _l_ n 

47.38 ± 0. 


.43 


43.37 


8 


15 


TVyf—l, rcn 

ivirK 5U 


Q 1 

bl 


0. 


.0234 






1.31 ± 


.23 


43. 16 


< 1 


4 


NGC 4355 


S2 


0. 


0084 


16.8 




i T o n 1 /"\ 

17.89 ± 0. 


31 


42.79 


40 


14 


M /■"! (~\ /lone 

NGC 4395 


Si. 8 


0. 


.0011 


3.6 




1.55 ± 0. 


.29 


40.38 


2 


18 


NGC 4507 


S2 


0. 


.0118 






10.93 ± 


.49 


43.48 


59 


15 


'M/"'/ - 1 A r no 

NGC 4593 


SI 


0. 


0090 


39.5 




5.86 ± 0. 


.25 


43.04 


< 1 


3 


NGC 4945 


S2 


0. 


0019 


5.2 




19.92 ± 


.36 


41.81 


220 


6 


LibO 323-G077 


Si. 2 


0. 


0150 






2.78 ± 0. 


.35 


43.09 


30 


8 


IbU J 13091 + 1137=JNGO 4992 




0. 


0251 




2 


3.49 ± 


.41 


43.65 


90 


13 


IGR J 13149+4422— Mrk 245 


S2 


0. 


0366 




7 


2.16 + 


.38 


43.77 






Cen A 


NLRG 


0. 


0018 


4.9 




56.09 ± 0. 


.32 


42.21 


11 


15 


MCG -6-30-15 


SI. 2 


0. 


0077 






3.62 ± 0. 


.38 


42.62 


< 1 


5 


Mrk 2 bo 


S2 


0. 


0399 






1.74 + 0. 


.30 


43.76 






ATT "1 O A A £> f\ 

4U 1344— bU 


SI 


0. 


0130 




6 


5.96 ± 


.25 


43.30 


< 1 


3 


t j^i ,1 o on A 

1C 4329A 


Si. 2 


0. 


.0160 






16.15 ± 


.51 


43.92 


< 1 


5 


Circinus galaxy 


S2 


0. 


0014 


4.2 




18.7b ± 0- 


26 


41.60 


400 


6 


NGC 550b 


Si. 9 


0. 


0062 


28.7 




13.34 ± 


66 


43. 12 


2.6 


15 


1GK J 14493— 5534=21VLAbX. J 14491283 — 553bl94 










8 


1.64 ± 0. 


.23 




10 


8 


ItrK J 14515 — 5542= W K.1S. 4<i74 


S2 


0. 


0180 




3 


1.55 ± 


.22 


43.00 


< 1 


8 


ICjH J 14552 — 5133— WKK 443o 


NLS1 


0. 


0160 




3,9 


1.40 ± 


.23 


42.85 


< 1 


4 


IC 4518 


S2 


0. 


0157 






2.44 + 0- 


.23 


43.08 






WKK 0092 


SI 


0. 


0156 






1.66 ± 0. 


.24 


42.91 


< 1 


19 


ICjH J lblo5 — 5925— WKK b471 


NLS1 


0. 


0350 




3,9 


1.73 ± 


.24 


43.64 






EbO 137-G34 


S2 


0. 


0092 






1.68 ± 0. 


.23 


42.45 






1GH J lb452 — 303b— 2MAbA. J lb4o!523 — 3035037 


SI 


0. 


0313 




6 


2.61 i 


.24 


43.72 


< 1 


8 


IN UL DZZl 




0. 


0050 


in 

19.4 




1.90 ± 0- 


.28 


41 .94 


1 


20 


1GK J10555 — 5203 


Si. 2 


0. 


0540 




3 


2.93 ± 


.21 


44.26 


< 1 


8 


Mrin conn 

IN GO b300 


S2 


0. 


.0037 


14.3 




4.71 ± 0. 


.45 


42.07 


25 


15 


T 1 T3 T 1 TO ("1 /I QEEyl 

lljil J17JU4 — J554 










10 


1.13 ± 


.17 




13 


2 1 


OrJrtb 1734 — 292 


SI 


0. 


.0214 






7.42 ± 0- 


. 14 


43.83 


< 1 


22 


Tr' O T17/I10 i o i o oTJ 1 i Ton i inin 

1GH J 17415— 1212— 2K 1739.1 — 1210 


SI 


0. 


0372 




11 


2.55 ± 


.29 


43.86 


< 1 


23 


1GK J174oo — o25o 


SI 


0. 


0200 




3 


1.91 ± 0. 


. 14 


43.19 


< 1 


8 


Idrt J1751J — 2U11 


Si. 9 


0. 


0470 




3 


2.32 ± 


.17 


44.03 






1GK JloU27 — 1455 


SI 


0. 


0350 




12 


2.93 i 


.22 


43.87 


< 1 


24 


qp onn 


DT T) /~1 


0. 


0561 






6.16 ± 0. 


.64 


44.61 


< 1 


3 


IGR J18559+1535=2E 1853.7+1534 


SI 


0. 


0838 




13,14 


2.27 ± 0. 


.23 


44.54 


< 1 


23 


1H 1934-063 


SI 


0. 


0106 






1.77 + 0. 


.31 


42.59 


< 1 


4 


NGC 6814 


SI. 5 


0. 


0052 


22.8 




4.73 ± 0. 


.41 


42.47 


< 1 


5 


Cygnus A 


NLRG 





0561 






5.77 + 0. 


.33 


44.58 


20 


26 


IGR J2018+4043=2MASX J20183871+4041003 










15 


1.89 + 0. 


.28 




7 


8 


IGR J20286+2544^MCG +04-48-002 


S2 


0. 


0142 




16 


3.31 ± 0. 


.58 


43.12 


50 


8 


Mrk 509 


SI. 2 


0. 


.0344 






5.51 ± 


.83 


44.13 


< 1 


15 


IGR J21178-|-5139=2MASX J21175311-|-5139034 










17 


1.93 ± 0. 


.34 








IGR J21247+5058 


SI 


0. 


0200 




12 


8.56 ± 0. 


.34 


43.84 


< 1 


24 


IGR J21277+5656 


SI 


0. 


.0144 




13 


2.68 ± 0. 


.45 


43.04 


< 1 


8 


NGC 7172 


S2 


0. 


0087 


33.9 




5.99 + 


.50 


42.92 


13 


15 


MR 2251-178 


SI 


0. 


0640 






4.75 + 0. 


.48 


44.62 


< 1 


27 


NGC 7469 


SI. 2 


0. 


0163 






4.74 + 0. 


.78 


43.40 


< 1 


3 


Mrk 926 


SI. 5 


0. 


0469 






3.56 + 


.53 


44.21 


< 1 


3 
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Table 1. -continued 



Object 




Class 3 




z 


D Notc b 




^17-60 kc 


V 


log -Lh x 




Ref c 


(sources at \b\ < 5° are given in 


bold) 








Mpc 


10 


erg/s/cm 


crg/s 


10 22 cm~ 2 














Blazars 














87GB 003300.9+593328 




Bl 


0.0860 






1.03 ± 0.15 


44.22 


< 1 


1 


S5 0836+71 




Bl 


2.1720 






3.47 ± 0.33 


48.03 


< 1 


1 


PKS 1219+04 




Bl 


0.9650 






1.30 ± 0.23 


46.74 


< 1 


4 


3C 273 




Bl 


0.1583 






13.83 ± 0.23 


45.92 


< 1 


1 


3C 279 




Bl 


0.5362 






1.50 ± 0.28 


46.18 


< 1 


1 


Mrk 501 




Bl 


0.0337 






4.53 ± 0.43 


44.02 


< 1 


1 


PKS 1830-211 




Bl 


2 


.5070 






3.35 ± 0.22 


48.17 


< 1 


25 


BL Lac 




Bl 


0.0686 






1.88 ± 0.38 


44.28 


< 1 


1 


3C 454.3 




Bl 


0.8590 






13.92 ± 0.55 


47.64 


< 1 


1 


Transiently detected AGN 


NGC 0526A 




SI. 5 





.0191 






3.72 ± 


0.97 


43.43 






ESO 297-G018 




S2 





.0252 






4.41 ± 


0.94 


43.75 






NGC 1052 




S2 





.0050 


17.8 




2.12 ± 


0.43 


41.91 






IGR J03334+3718 




SI 





.0547 


18 




1.96 ± 


0.41 


44.09 






LEDA 168563 




SI 





.0290 






4.42 i 


1.11 


43.88 






AKN 120 




SI 





.0327 






8.38 ± 


3.52 


44.26 






MCG 8-11-11 




SI. 5 





.0205 






5.61 ± 


1.51 


43.67 






IRAS 05589+2828 




SI 





.0330 






3.03 ± 


0.62 


43.83 






LEDA 096373 




S2 





.0294 






2.71 ± 


0.61 


43.68 






PG 0804+761 




SI 





.1000 






1.21 ± 


0.33 


44.43 






NGC 4051 




SI. 5 





.0023 


17.0 




2.11 i 


0.54 


41.87 






NGC 4138 




SI. 9 


o 


.0030 


17.0 




2.39 ± 


0.49 


41.92 






WAS 49B 




S2 





.0630 






0.95 ± 


0.35 


43.91 






NGC 4253 




SI. 5 





.0129 






1.53 ± 


0.31 


42.70 






NGC 4258 




SI. 9 





.0015 


6.8 




1.90 ± 


0.56 


41.03 






XSS J12389-1614=IGR J12391 


-1612 


S2 





.0367 


19 




3.01 ± 


0.62 


43.92 






WKK 1263 




S2 





.0244 






1.28 ± 


0.31 


43.19 






PKS 1241-399 




QSO 





.1910 






1.77 ± 


0.40 


45.21 






Mrk 783 




SI. 5 





.0672 






1.23 ± 


0.32 


44.08 






IGR J13038+5348=MCG 09-21-096 


SI 





.0299 


18 




2.11 ± 


0.45 


43.58 






NGC 5033 




SI. 9 





.0029 


18.7 




1.21 ± 


0.28 


41.71 






ESO 383-G018 




S2 





.0124 






1.74 ± 


0.38 


42.72 






NGC 5252 




SI. 9 





.0230 






2.60 ± 


0.78 


43.44 






IGR J14175-4641 




S2 





.0760 


3 




1.30 ± 


0.28 


44.21 






NGC 5548 




SI. 5 





.0172 






2.06 ± 


0.50 


43.08 






ESO 511-G030 




SI 





.0224 






3.33 ± 


0.91 


43.53 






IGR J14471-6319 




S2 





.0380 


3 




0.95 ± 


0.23 


43.45 






NGC 6240 




S2 





.0245 






4.66 ± 


1.38 


43.75 






3C 382 




SI 





.0579 






3.96 ± 


1.75 


44.45 






ESO 103-G035 




S2 





.0133 






6.64 ± 


1.34 


43.37 






XSS J19459+4508=IGR J19473- 


f4452 


S2 





.0532 


20 




1.52 ± 


0.45 


43.96 






3C 403 




NLRG 





.0590 






0.96 ± 


0.33 


43.85 






4C +74.26 




BLRG 





.1040 






3.93 ± 


1.22 


44.98 






S5 2116+81 




SI 





.0840 






2.91 ± 


0.98 


44.65 






NGC 7314 




SI. 9 





.0048 






2.17 ± 


0.58 


41.99 






Mrk 915 




SI 





.0241 






2.62 ± 


0.56 


43.49 







a Optical AGN class: Bl - blasar (BL Lac object or flat-spectrum radio quasar), SI, SI. 2, SI. 5, SI. 8, SI. 9, S2 — Scyfcrt galaxy, NLS1 — narrow-line 
Seyfert 1 galaxy, BLRG - broad-line radio galaxy, NLRG — narrow-line radio galaxy, QSO — quasar. 

b Note or reference for INTEGRAL discovered AGN: (1) IBurenin eTal] i2006bl) . (2) ISazonov e t alj I200 5|). (3 1 IMasetti et"ai] I2006dl) . (4) 
galaxy=lRXS J0944 36. 5-263353, a Seyfert 1 nucleus is suggested by the 6dF spectrum. (51 IMorelliTt~aiTT200a) . =SWIFT J1038. 8-4942 (6) 

|Mase^ti_et_alJ |2006j), (7) SPSS s pectru m indicates a Seyfe rt 2 n ucleus, (8) Swi f t loca lizatio n consistent with this galaxv^(jj1 jReyrii\^evetalJ 
120091. (101 IBassani et alj I200 5y. fill ITorres et all 1200411. (12) IMasetti et al J 1200411. (131 IBikmaev et alj I2006D . (14) IMasetti eH!dTT27)06bll 
151 Ikenn ea et alJ 1 20031. (161 lMase.tti^r^jTl2T)06cl) . (171 IBassani et alj f 20o 3l . (Ts l IBureninTt^jTIT oOfill . (191 =2MASX J12390630-161 0472 
Sazonov et^dl2005t IMasetti eta!l200r^ . (201 =2MASX J19471 938+4449425 jSazonov et^lTl2ToalBikmaev et al.120'071 IMasetti et al.l2006bl) . 
c Quoted Af H v alue o r li mit is adopted from or based on: (l)IDonatoet all l2005l) , (2) RXTE data, (3) ASC A data, (4) ROSAT data, (5) 

[ Reynolds et alj 1199711 . (61 IMatt et al. I 1200011 . (71 IMalizia et alj Il997l1 . (81 Swift data, (91 IRisaliti et alj 1200511. Nm is str ongly variable, (10) 

ILewrsTt"al.l 1200511 "(111 IVignali et alj 119931 . < 121 ILmnlerT^T] 1200311 complex X-ray absorption. ( 131 ISazonov et at] 1200511 . ( 141 IRisaliti et all 
|200^)", (15) ^azono^^R^ni^s^3'"l^^^)7( 16l Swift data, it is also possible that the X- ray sp ectrum is reflection-dominated with A^h > 10 
l Aiello et alJl20061). (17) | V i gnaliTtafl " ^2002ll . (18l lMoran et alj 120051) . (19l lWoudt et alj 11993) . (20) ILevenson et alj 120011). (21) IBassani et all 
1200511. (22) Sazonov et all feo04bl). (231 E i nstein data, (24) assumed upper limit based on Seyfert 1 classification, (25) IQshima et^ilTl20UlTl . (26) 

lYoung et alj 1200211 . (27) IReeves fc Turner! 1200011 . 



We also collected information on the X-ray absorption 
columns for the AGN of the main sample (column 8 of 
Table For this purpose we either used the literature 
or analyzed publicly available spectral data from different 
X-ray instruments, including ASCA/GIS, Einstein/IPC, 
RXTE/PCA, and Swift /XRT. The sp ectral modelling was 
done with XSPEC vl2 l|Arnaudlll99(| . The X-ray spectra 
(usually above 2 keV) were fitted by a model consisting 
of a power law photoabsorbed by a column TVh of neu- 
tral material and a fluorescent neutral-iron emission line 



at rest energy 6.4 keV if the latter was required by the 
fit (model zphabs x (powerlaw + Gaussian) in XSPEC). If 
the quality of the data allowed, we considered the power- 
law slope a free parameter; other wise, the photon in- 
dex was fixed at th e canonical (e.g. iRevnolds et al.lll997t 
iTurner et, all Il997fl AGN value T = 1.8. Note that we 
did not try to evaluate Nn if it was clear that iVn < 
10 cm -2 (after subtraction of the Galactic absorption), 
as we regard such sources as unobscured throughout the 
paper. Finally, for several Seyfert 1 galaxies in our sam- 
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pie, which are expected to have little X-ray absorption, we 
used ROSAT/PSPC soft X-ray data to verify that indeed 
A^h < 10 22 cm -2 for them. 

2.1. Unidentified sources 

There are 24 unidentified sourc es among those detec ted 
on the average IBIS /ISGRI map ( Kr ivonos et al . 2006), so 
our all-sky AGN sample can be significantly incomplete. 
However, since most of these objects are located close to 
the Galactic plane, we can optimize our statistical analysis 
by excluding the Galactic plane region \b\ < 5° from con- 
sideration. This leads to only a small decrease in the num- 
ber of identified AGN - to 74, but to a dramatic reduction 
of the number of unidentified sources - to 7. The coordi- 
nates and fluxes of these high Galactic-latitude sources - 
AGN candidates - are given in Tabled We note that the 
apparently high sky density of unidentified sources in the 
direction I ~ 330°, b ~ 10° is probably the result of the 
extensive INTEGRAL observations of the Galactic plane 
and of such interesting sources as Cen A, NGC 4945, and 
SN 1006. 

The treatment that follows is thus based on the |6| > 5° 
AGN sample. Furthermore, we restrict our consideration 
to the emission-line AGN (66 in total) and so do not dis- 
cuss blazars anymore. 

3. X-ray absorption distribution 

We first address the X-ray absorption distribution of 
nearby (z <; 0.1) emission-line AGN. One of the 
most important findings of the RXTE Slew Survey 
was that the fraction of significantly X-ray absorbed 
AGN (log Ah > 22) drops with increasing luminosity 
I Sazonov fc R.evnivtsevlEoO^ . Motivated by this result, 
we divided our current INTEGRAL sample into two parts: 
AGN with logihx < 43.6 and those with logLhx > 43.6, 
where L nx is the luminosity in the 17-60 keV band. The 
chosen dividing luminosity approximately corresponds to 
the value of the 3-20 keV luminosity (logL x = 43.5) 
that we used in the RXTE case given that typically 
log(Lh x /i x ) ~ 0.1 for local AGN, as will be shown in H4.3I 
below. This boundary also approximately corresponds to 
the bend of the hard X-ray luminosity function (@. The 
resulting low- and high-luminosity samples include 42 and 
24 objects, respectively. 

Figure shows the observed Ah distributions for these 
subsamples. It can be seen that, while ~ 66% of local 
logLhx < 43.6 AGN are obscured (logA^ > 22), the cor- 
responding fraction is only ~ 24% among the higher lu- 
minosity AGN. This result confirms our XSS finding but 
it is obtained here in a more straightforward manner. In 
the INTEGRAL case, the efficiency of source detection 
at energies above ~ 20 keV is practically independent of 
intrinsic source absorption, as long as logA^H < 24.5. In 
the XSS case, the detection was based on 3-20 keV fluxes, 
so we had to correct the observed Ah distribution for the 
significant loss in sensitivity to sources with log Ah > 23. 



15 - 



10 



log L<43.6 



3 
55 



15 



10 - 



log L>43.6 



22-22. 5 23-23.5 24-24.5 unknown 

<22 22.5-23 23.5-24 >24.5 



Log N H (cm- a ) 



Fig. 2. Observed X-ray absorption distribution of the low- 
luminosity AGN (top panel), and high-luminosity AGN 
(bottom panel). The shaded part of each diagram shows 
the number of AGN with unknown A%. 

We note that the above conclusion with respect to the 
fraction of obscured AGN as a function of luminosity is 
unlikely to change once the absorption columns are mea- 
sured for those INTEGRAL AGN for which this informa- 
tion is not available yet, since there are only 4 and 3 such 
sources out of the 42 low- and 24 high-luminosity AGN, 
respectively (see Fig- - 

We point out that the observed (see Fig. EJ small frac- 
tion of substantially Compton-thick AGN - sources with 
log Ah ^ 24.5 - should be interpreted with caution, since 
in their case the intrinsic hard X-ray flux is expected to be 
substantially suppressed by Comptonization on the cold 
material obscuring the AGN. The discussion of the frac- 
tion of Compton-thick AGN in the local Universe will be 
continued in Sect. |SJ 

4. Hard X-ray luminosity function 

We now define the hard X-ray luminosity function of 
nearby emission-line AGN as 4>{L\^) = dN agn / 'dlog £hx- 
We first estimated d>(L hv) in binned form using the 1/V m 
method i.e. by summing l/V m (Lhx,i) val- 

ues for our AGN in specified luminosity intervals, where 
Vm(Lhx,i) is the space volume in which an AGN with lumi- 
nosity Li lx .i could be detected by the survey. The V m calcu- 
lation was based on the dependence of the sky area covered 
by the survey on the achieved sensitivity, which is shown in 
Fig. We note that - 75% and ~ 50% of the \b\ > 5° sky 
have been covered down to 5 and 3 mCrab, respectively, 
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Table 2. Unidentified INTEGRAL sources at \b\ > 5° 



Source 


RA 


Dec. 


1 


b 


-^17-60 kcV 


Likely associations 




(J2000.0) 






10 -11 crg/s/cm 2 




IGR J02466-4222 


41.65 


-42.37 


253.48 


-62.07 


3.08 ± 0.55 




IGR J09522-6231 


148.05 


-62.52 


283.83 


-6.50 


1.17 ± 0.21 




IGR J13107-5551 


197.68 


-55.87 


305.66 


6.90 


1.85 ± 0.30 


XMMSL1 J131042. 6-555206 (Read et al. 20051=PMN J1310-5552 


IGR J14561-3738 


224.04 


-37.65 


328.61 


18.93 


1.40 ± 0.26 


XSS J14562-37 


IGR J16500-3307 


252.51 


-33.10 


349.73 


7.34 


1.61 ± 0.23 


1RXS J164955. 1-330713 


IGR J16562-3301 


254.06 


-33.03 


350.60 


6.37 


1.97 ± 0.21 


SWIFT J1656. 3-3302 fTueller et al. 20061 


IGR J17350-2045 


263.75 


-20.76 


5.69 


6.36 


1.29 ± 0.18 






10- 11 10- 10 10- 9 40 41 42 43 44 45 

Limiting flux (17-60 keV), erg s _1 cm -2 Log 17-60 keV luminosity (erg s _1 ) 



Fig. 3. Fraction of the total sky area as a function of sur- 
vey sensitivity (the dashed line). The solid line shows the 
corresponding fraction for the \b\ > 5° sky 

where 1 mCrab corresponds to ~ 1.4 x 10~ n erg s cm -2 
(17-60 keV). The corresponding fractions are somewhat 
higher for the whole sky due to the overexposure of the 
Galactic plane region by INTEGRAL. 

We performed the calculation in the luminosity range 
40 < logihx < 45.5, which includes all of our emission- line 
AGN, 66 objects in total. The resulting binned luminosity 
function with estimated la statistical errors is shown in 
Fig.H 

We next sought an analytic representation of the lumi- 
nosity function in the form of a smoothly broken power- 
law model, as is conventional in AGN research: 



A 



(Lhx/£*) 71 + (W-M 72 



(1) 



For this purpose we used the maximum likelihood estima- 
tor 



L = -2 ^ In 



</ , (^hx,i)V r m (ihxa) 
/ 0(ihx)l / m(ihx)dlog£hx' 



(2) 



Fig. 4. Hard X-ray luminosity function of local emission- 
line AGN obtained with INTEGRAL in binned form (solid 
circles and la error bars). The solid line shows the analytic 
approximation given by Eq. and Table 



where i runs over all sampled AGN. 

By minimizing L, we found the best-fit values for the 
break luminosity L*, as well as for the slopes 71 and 72. 
Since the maximum likelihood method directly determines 
only the shape of the function, we derived the normaliza- 
tion A by requiring that the number of AGN predicted by 
the model be equal to the number of AGN in the sam- 
ple. The resulting best-fit model is shown in Fig. 0] and 
its parameters are given in Table together with their 
estimated la errors; the error for the coefficient A is not 
quoted, since this parameter is strongly correlated with 
the others. The analytic fit is apparently in good agree- 
ment with the binned 0(L n x), and the good quality of the 
fit is confirmed by the Kolmogorov-Smirnov test. 
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Table 3. Parameters of the hard X-ray luminosity func- 
tion 



Parameter 


Value and la range 


30 




logL, 






43.40 (43.12 — 43.68) 






7i 






0.76 (0.56-^0.94) 






72 






2.28 (2.06-^2.56) 


20 




A (Mpc -3 ) 






3.55 x 10" 5 






niT-eo kov(> 40) (1(T 3 


Mpc" 


- 3 ) 


9 (4 -=-18) 


10 




nn-ao k c v(> 41) (HT 3 


Mpc" 




1.4 (0.9^2.0) 






/917-60 kcv(> 40) (10 38 


erg s~ 


1 Mpc -3 ) 


14.1 (11.8 -r 17.1) 







P17-6O kov(> 41) (10 38 


erg s~ 


1 Mpc" 3 ) 


12.4 (11.0 ~ 14.0) 




Pks 


> 0.9 







1 - 



4.1. Cumulative number density and luminosity density 

By integrating the derived luminosity function over lumi- 
nosity, we can find the total number density of nearby 
emission-line AGN with log Lhx > 40: 



n(> 40) 



c/>(L hx ) d log L hx - 



(3) 



40 



The resulting value is given in Tabled together with its la 
uncertainty range. It should be noted that this cumulative 
number density is strongly dominated by AGN at the low 
end of the probed luminosity range and that there is only 
one AGN in our sample (NGC 4395) with logL hx < 41. 
For this reason we also provide our estimate of (similarly 
defined) n{> 41) in TableEl 

We can also determine the cumulative luminosity den- 
sity of nearby emission-line AGN with log Lhx > 40: 



P17-6O kcv(> 40) 



4>(L h ^)L^dlogL 



hx • 



(4) 



The result is presented in Tabled Also given in the table is 
the value of pi7_eo kev(> 41), which is better constrained 
by the data. 

4.2. Systematic uncertainties 

There are some systematic effects that may affect the lu- 
minosity function obtained. First, our input AGN sam- 
ple may be incomplete because there are 7 unidentified 
INTEGRAL sources at \b\ > 5°, some of which may even- 
tually prove to be emisison-line AGN. Assuming that there 
is no strong luminosity bias for these sources, the ampli- 
tude of our luminosity function can be underestimated due 
to incompleteness by ~ 10% at most. 

Another cause of concern is that 15 of our \b\ > 5° 
emission-line AGN have been at least once the target of a 
pointed INTEGRAL observation, i.e. we are not dealing 
with a truly serendipitous survey. In order to assess the 
maximum possible systematic effect associated with these 
AGN pointings, the following argument can be suggested. 
Suppose that an area S of the sky was covered by the 
survey down to flux f\ or lower, while the total area of 
the survey is So- Suppose next that a total of N sources 
with flux /1 < / < / 2 were detected, of which N t were 



s 



0.9 



o.e 



10- 1 



10- 10 

Flux (17-60 keV), 



10- 9 



erg s 



Fig. 5. Top: Differential distribution of hard X-ray fluxes 
for all INTEGRAL emission-line AGN at \b\ > 5° (upper 
histogram) and only for those targeted by INTEGRAL 
(filled histogram). Bottom: The ratio of the minimum and 
maximum estimates of the sky density of sources accord- 
ing to Eq. JSJ. 



observational targets. Then the sky density p of sources 
with fi < / < f<i can be expected to be bound in the 
range (p min < p < p ma x): 



N - N t N t N 

c h o~ < P < "o- 



(5) 



Our calculation of the luminosity function was im- 
plicitly based on using p max . Figure [S] shows the ratio 
Pmin/Pmax as a function of flux for our AGN survey. It 
can be seen that the effect under consideration is strongest 
at medium fluxes (~ 2-6 10 -11 erg s _1 cm~ 2 ) where the 
fraction of targeted AGN is relatively high, while the sky 
coverage is relatively low. However, even at the maxi- 
mum the effect is less than 17%. This value can be con- 
sidered an upper limit on the possible overestimation of 
the luminosity function due to the presence of targeted 
AGN in the sample. We note that implementing observa- 
tions of "empty" extragalactic regions with INTEGRAL 
(Proposal 0320108) has allowed us to significantly reduce 
the non-randomness of the survey. 

We may conclude from the above that the total system- 
atic uncertainty in the luminosity function due to source 
identification incompleteness and AGN pointings is proba- 
bly <; 10% (taking into account that the two effects tend to 
counterbalance each other), which is less than the present 
statistical uncertainties. 
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Fig. 6. V/V m ratio averaged over luminosity bins for the 
AGN sample used for the construction of the luminosity 
function. The error bars represent la statistical uncertain- 
ties. 



Finally, the standard V/V m test l)Schmiddll968|) ap- 
plied to our AGN sample demonstrates (Fig. HJ) that the 
spatial distribution of the studied AGN population is con- 
sistent with being homogeneous. 

4.3. Comparison with the 3-20 keV luminosity function 

It is interesting to compare the local AGN luminosity func- 
tion obtained with INTEGRAL in the hard X-ray band 
(17-60 keV) with that previously constructed in the ad- 
jacent, medi um X-ray band (3-20 keV) u sing the RXTE 
Slew Survey (|Sazonov fc RevnivtsevlEooil) . 

For this comparison we need to know how the 17- 
60 keV luminosities (Lhx) of Seyfert galaxies are related to 
their 3-20 keV luminosities (L x ). One possibility is to di- 
rectly compare L^ x with L x for those emission-line AGN 
detected during both the XSS and the INTEGRAL all- 
sky survey. We have in total 25 such cases, for which we 
show in Fig.[7|L x vs. L^ x and the ratio L x /L^ x as a func- 
tion of Ah- These data are compared with the expected 
dependence of L x /Lhx on Ah assuming that the intrinsic 
(unabsorbed) AGN spectrum consists of a power law with 
a photon index L = 1.8 plus a Compton reflection compo- 
nent (XSPEC model pexrav. lMaedziarz fc Zdziarskll995[) 
of different relative amplitudes R (= f2/27r, for a zero incli- 
nation angle) . As expected (see the upper panel of Fig. EJ , 
for most of the AGN L x < L^ x mainly due to the pres- 
ence of intrinsic X-ray absorption and (to a lesser degree) 
of Compton reflection. The measured L x /Lhx ratios are 
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Fig. 7. Top: The 3-20 keV luminosity as a function of 
17-60 keV luminosity for AGN detected during both the 
RXTE Slew Survey and INTEGRAL all-sky survey. The 
line indicates equality of the two luminosities. Bottom: 
The 3-20 keV/ 17-60 keV luminosity ratio as a function 
of Nji (points). The data shown at logA^H = 21 represent 
unobscured AGN (log Ah < 22). Also shown are expected 
luminosity ratios in the absence of source variability for 
an intrinsic (unabsorbed) spectrum consisting of a power 
law (r = 1.8) and a Compton reflection component with 
relative amplitude R — (dotted line), R — 0.5 (dashed 
line), and R — 1 (solid line). 



apparently consistent with the expected trend, some- 
what more so in the case of a non-negligible reflection 
component (R ~ 0.5-1). The observed significant (a fac- 
tor of ~ 2) scatter of the data around this trend can be 
naturally explained by intrinsic variability of the sources, 
since the RXTE and INTEGRAL observations are sepa- 
rated by years. 

Taking the outcome of this test into account and con- 
sidering the abundant literature on the broad-band X- 
ray spectra of Sey fert galaxies (e.g. iGondek et aO Il996t 
IPerola et alJl2002|) . we may adopt to a first approxima- 
tion that local emission-line AGN have a universal in- 
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Fig. 8. The best-fit model of the hard X-ray AGN lumi- 
nosity function obtained with INTEGRAL (line) in com- 
parison with the luminosity functi on based on the RXTE 
Slew Survey in the 3-20 keV band ( Sazonov & Revnivtsevl 
2004), here converted to the 17-60 keV band (points with 
error bars). 



trinsic spectrum consisting of a power law with T = 1.8 
and a substantial reflection component (R <~ 0.5-1). We 
can then estimate the average ratio (L x /Li lx ) for low- 
(logLhx < 43.6) and high- (logLhx > 43.6) luminosity 
AGN taking their (different) observed A% distributions 
(Fig. EJ) into account. For the low-luminosity AGN, we 
find (L x /L hx ) = 0.79 (0.67) for R = 0.5 (1). For the 
high-luminosity AGN, most of which are unobscured, the 
corresponding values are 0.96 (0.84). Note that, in this 
context, changing the reflection amplitude R by 0.5 corre- 
sponds to changing the power-law slope by ~ 0.1. Given 
the above numbers, we may adopt to a first approxima- 
tion that log(Lhx/-^x) ~ 0.1 regardless of luminosity, and 
accordingly convert the XSS luminosity function to the 
17-60 keV band. 

The result is shown in Fig. [S] in comparison with 
the hard X-ray luminosity function measured directly by 
INTEGRAL. Note that the amplitude and error bars for 
the XSS luminosity function take into account the sys- 
tematic uncertainty d ue to possible incompleteness of the 
XSS sample of AGN llSazonov & Revnivtsevl Eofil . The 
two luminosity functions agree with each other, as is for- 
mally confirmed by a comparison of the parameters of 
the best-fit analytic models. However, there is an indi- 
cation that the XSS luminosity function lies somewhat 
below the INTEGRAL one. Here it is best to consider 



the cumulative luminosity density. With RXTE we mea- 
sured P3-20 kcv(> 41) = (5.2 ±1.2) x 10 38 erg s^ 1 Mpc" 3 , 
which for the {L\ lyi / L x ) ratio adopted above implies that 
P17-60 kov(> 41.1) = (6.5 ± 1.5) x 10 38 erg s" 1 Mpc" 3 . 
With INTEGRAL we found p l7 _ m kcV (> 41) = (12.4 ± 
1.5) x 10 38 erg s _1 Mpc" 3 , i.e. a factor of 1.9 ±0.5 higher 
luminosity density. 

5. Implications for the cosmic X-ray background 

It has become a common paradigm that the bulk of the 
cosmic X-ray background (CXB) is composed of emis- 
sion from all AGN in the Universe. This conclusion is 
mainly based on the fact that deep extragalactic X-ray 
surveys h ave resolved ~ 80% of the C XB at energies below 
~ 8 keV (|Hickox fc Markevitchll2006h . However, the max- 
imum of the CXB vF v spectrum is located at ~ 30 keV, 
and at these energies not more than ~ 3% of the CXB has 
been resolved into point sources in the deepest exposures 
with INTEGRAL (and much less by previous missions). 
Therefore, resolving the CXB near its peak will be one 
of the main tasks of future hard X-ray missions. For the 
time being, it is interesting to address the following ques- 
tion: is the X-ray absorption distribution measured with 
INTEGRAL for the local AGN population consistent with 
that required for distant quasars to explain the CXB spec- 
trum? 

To answer this question, we can deduce relative frac- 
tions of AGN with different absorption columns from the 
measured local A% distribution (Fig.p2J|: 



logA^H < 22.0 

22.0 < logA^H < 22.5 

22.5 < logA^H < 23.0 

23.0 < logA^H < 23.5 

23.5 < logA^H < 24.0 

24.0 < logA^H < 24.5 

logA^H > 24.5 



0.37 ±0.10, 
0.20 ±0.07, 
0.05 ±0.04, 
0.13 ±0.06, 
0.16 ±0.07, 
0.07 ±0.05, 
0.02 ±0.02. 



(6) 



These fractions take into account that low- (logLi lx < 
43.6) and high- (logLi lx > 43.6) luminosity AGN con- 
tribute ~ 90% and ~ 10%, respectively, to the local lumi- 
nosity density. The quoted errors result from Poissonian 
statistics and we have ignored the systematic uncertainty 
associated with the lack of iVn measurements for several 
AGN in our sample. 

We next adopt a fiducial intrinsic AGN spectrum con- 
sisting of a cutoff power-law component and the corre- 
sponding Compton reflection component with relative am- 
plitude R: 



dN 7 
~dE 



AET r exp( 



-E/Ef) + Rf(E), 



(7) 



where f(E) is again descri bed by the pexrav model 
l|Magdziarz fc Zdziarskil ll995A. We adopt L = 1.8, E { = 
200 keV, and R = 0.5 or R = 1 as reference values. 

We can now build the composite spectrum of the lo- 
cal AGN by propagating the intrinsic spectrum given by 



10 



Sazonov et al.: Hard X-ray luminosity function of nearby AGN 



Eq. (JJJ through different absorption columns and sum- 
ming up the resulting spectra with the weights introduced 
in Eq. ©. The resulting spectrum is shown in Fig. El Note 
that we have neglected the contribution of substantially 
Compton-thick AGN (logA^H > 24.5) since it is expected 
to be small (~ 3%) given the observed JVh distribution 
(Fig. EJ , whereas the spectra of Compton-thick AGN can 
assume quite different shapes depending on t he geometry 
of th e obscuring and reflecting material (e.g. iMatt et alJ 
1200(1) . 

The composite spectrum shown in Fig. [5] has been 
normalized so as to reproduce the local AGN luminos- 
ity density measured with INTEGRAL, pnso kcv(> 40) 
(Table El, i.e. this spectrum represents the volume emis- 
sivity of all local emission- line AGN with logLhx > 40. In 
Fig. E| we indicate the range of uncertainty in the compos- 
ite spectrum resulting from the uncertainties in the A?h 
distribution and in the pn~eo kcv(> 40) value. It can also 
be seen that the difference between the normalized com- 
posite spectra for R = 0.5 and R = 1 is only noticeable 
below ~ 15 keV. 

Suppose now that the same spectral shape character- 
izes AGN at all redshifts and that only the AGN lumi- 
nosity density evolves with z. Then the cumulative AGN 
emission obser ved at z = will have a spectrum (for a flat 
cosmology, e.g. lSazonov et alll2004|) 



1(E) = 



4nHn 



e(z)F((l+z)E) 



oJo (i + z )[n m (i + z) 3 + n A ] 



1/2 



dz, 



(8) 



where F(E) is the template spectrum shown in Fig.[5]and 
e(z) is a function describing the evolution of the luminosity 
density. 

The best constraints on e(z) have so far been 
coming from deep X-ray surveys with Chandra and 
XMM-Newton. In Fig. we reproduce the results of 
iBarger et al.l l(2005h on the redshift dependence of the rest- 
frame 2-8 keV luminosity density of AGN, which reflects 
the observed pure luminosity evolution of the AGN lumi- 
nosity function at z <J 1.5. Between z ~ 0.2 and z ~ 1 
the evolution is well-constrained and is consistent with a 
power-law dependence (z a ) with a slope a ~ 3.2. At z ;> 1 
there is still significant uncertainty due to the incomplete- 
ness of source identification and redshift determination, 
and the evolution function e(z) is bound by the two sets 
of data points in Fig. 1101 or approximately by the following 
two limiting functions, also shown in the figure: 



e\(z) oc 



e 2 (z) 



(1 + z) 3 - 2 , 
ei (l)/z, 

(1 + z) 3 - 2 , 
e 2 (l), 



z < 1 
z > 1 

z < 1 
z > 1. 



(9) 



(10) 



By convolving our locally-determined composite AGN 
spectrum (Fig. |5J) with these two alternative evolution 
laws, we obtain the cumulative AGN spectra shown in 
Fig. El The shaded regions around these spectra repre- 
sent the superposition of the present uncertainties in the 



2 




0.01 



10 100 

Energy, keV 



1000 



Fig. 9. Composite spectrum of the local emission-line 
AGN obtained assuming that the intrinsic AGN spectrum 
is given by Eq. with r = 1.8, E { = 200 keV, and 
R = 0.5 (thick solid line), and taking the observed Nn dis- 
tribution of nearby AGN (Fig. EJinto account. The spec- 
trum is normalized to the measured local AGN 17-60 keV 
luminosity density. The shaded area indicates the uncer- 
tainty in the composite spectrum due to the uncertainties 
in the P17-60 kcv(> 40) value and in the Afe distribution. 
Also shown are the contributions to the composite spec- 
trum of AGN with different degrees of obscuration (curves 
with labels indicating logA^)- The dashed line shows the 
composite spectrum in the case of R = 1 . 

local P17-6O kcv(> 40) value and in the local Ah distri- 
bution. For comparison in Fig. ^2 the CXB broad-band 
spectrum is shown as recently measured with the d iffcr- 
ent instruments on INTEGRAL l|Churazov et al.l2006j) . It 
can be seen that the predicted cumulative AGN spectrum 
agrees very well, both in amplitude and in shape, with 
the CXB for both limiting scenarios of AGN evolution at 
z > 1. It should be noted, though, that our calculation 
was based on the INTEGRAL estimate of the cumulative 
luminosity density of local AGN with logihx > 40, so 
the inclusion of even lower luminosity AGN and normal 
galaxies could somewhat increase the normalization of the 
predicted cumulative AGN spectrum. 

We conclude that the combination of data from the 
INTEGRAL all-sky survey and deep extragalatic X-ray 
surveys seems to be able to naturally explain the ob- 
served CXB spectrum. Furthermore these data are con- 
sistent with the cosmological evolution of AGN at z <J 1.5 
(where the bulk of the CXB is produced) having occurred 
mostly in luminosity and much less (if any) in spectral 
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Fig. 10. Evolution with redshift of the rest-frame 2-8 keV 
luminosity density b ased on Chandra e xtragalactic sur- 
veys (adopted from iBarger et alJ l2005h . The measure- 
ments shown by the filled squares take only spectroscop- 
ically identified AGN into account, while those shown by 
the empty squares also take the maximal incompleteness 
into account. We assume that the true evolution is bound 
between the two limiting functions given by Eqs. and 
ifTTij) . shown by the dashed and solid lines, respectively. 

properties such as the X-ray absorption distribution (for 
a given value of Lhx/ L*(z)). 



Fig. 11. Predicted CXB spectrum based on the local AGN 
composite spectrum (Fig. for R — 0.5) and the AGN 
redshift evolution measured with Chandra (Fig. I10|). The 
upper thick solid line corresponds to the scenario of flat 
AGN luminosity density evolution at z > 1, the hori- 
zontally shaded region indicating the corresponding un- 
certainty. The lower solid line and the vertically shaded 
region correspond to the scenario of oc 1/z evolution at 
z > 1. The other presented curves show the contributions 
of AGN with different degrees of obscuration (the labels 
indicate log Ah). The points with error bars show the CXB 
spectrum measured with the JE T-X, IBIS/ISGRI, an d 
SPI instruments on INTEGRAL l)Churazov et al.ll2006|) . 



6. Conclusions 

We have used the INTEGRAL all-sky hard X-ray sur- 
vey to study some key properties of the local (z <; 0.1) 
AGN population. Since our source detection was based on 
17-60 keV fluxes, the survey is equally sensitive to AGN 
with X-ray absorption columns up to several 10 24 cm~ 2 , 
i.e. well into the Compton-thick domain (Ah > 1.5 x 
10 24 cm- 2 ). 

One of the most surprising results of the survey is that 
very few Compton-thick AGN have been detected, and 
all of them were known before. The observed fraction of 
Compton-thick objects is only ~ 10%. This estimate may 
increase to at most ~ 20% once the currently missing ab- 
sorption columns are measured. Clearly Compton-thick 
AGN are fairly rare. A note of caution is necessary here 
concerning very Compton-thick objects (log Ah ~ 25-26). 
Since their observed hard X-ray luminosities are expected 
to be <J 10% of their intrinsic luminosities, such sources 
may constitute a major fraction of AGN at luminosities 



near our survey's effective limit (Lhx ~ 10 40 -10 41 erg s *) 
and below. 

Another major result of our study is confirmation of 
the result of the RXTE Slew Survey that the local AGN 
absorption distribution is very different at low (log Lhx < 
43.6) and high (logLhx > 43.6) luminosities. While the 
fraction of obscured (log Ah > 22) objects is ~ 70% 
among the low-luminosity AGN, it is only ~ 25% among 
the high-luminosity ones. We note that a similar result 
was recently reported based o n a combination of different 
HEAO-1 (2-10 keV) surveys l)Shinozakill2006|l and is also 
emerging from the all-sky hard X -ray survey carried out 
by Swift l|Markwardt et ai]l2005|) . Furthermo re, a similar 
trend has been found at higher redshifts (e.g. lUeda et alJ 
l2003h . 

We measured the hard X-ray luminosity distribu- 
tion of local emission-line AGN. Its broken power-law 
shape is in good agreement with recent determina- 
tions in softer energy bands, in particular at 3-20 keV 



12 



Sazonov et al.: Hard X-ray luminosity function of nearby AGN 



(ISazonov fe Revnivtsevll2004|) and at 2-10 keV ijShinozakil 
2006| . We also found the cumulative luminosity density of 
AGN with L hx > 10 41 erg s" 1 to be pi 7 -60 kcv(> 41) = 
(12.4 ± 1.5) x 10 38 erg s^ 1 Mpc~ 3 . This is a factor of - 1.9 
higher than the XSS measurement in the 3-20 keV band 
assuming that typically (£17-60 kcv/^3-20 kcv) ~ 1-25. 
Thus the two determinations differ by — 2.8c, although 
there is a significant uncertainty in the average ratio 
(Li7_ 6 o kcv/^3-20 kcv)- Our INTEGRAL estimate of the 
AGN cumulative luminosity density in the 17-60 keV 
band also seems to be con sistent with the 2-10 keV es- 
timate of IShinozakil l)2006j) : p 2 -io kcv(> 42) = (7.2 ± 
1.4) x 10 38 erg s" 1 Mpc" 1 (for H = 75 km s" 1 Mpc" 1 ), 
given that on average Li7_60 kcv/^2-io kcv ~ 2. Finally, 
the hard X-ray luminosity function obtained in this 
work a grees with the one constructed bv lBeckmann et all 
l|2006|) based on a smaller sample of INTEGRAL AGN, 
improving it in terms of accuracy and reducing the bias 
associated with AGN pointed observations. 

In summary, the all-sky X-ray and hard X-ray surveys 
performed with RXTE, INTEGRAL, Swift, and HEAO-1 
have provided an accurate census of nearby unobscured 
and obscured AGN. This new information is not only in- 
teresting in its own right but also provides a reliable z = 
point for studies of AGN cosmological evolution and the 
growth of massive black holes. 

As a first attempt to apply the local AGN statistics 
in a cosmological context, we demonstrated that the spec- 
tral shape and amplitude of the CXB allow for the pos- 
sibility that the absorption distribution of AGN (for a 
given Lh K /L*(z) ratio) has not changed significantly since 
z — 1.5, while the AGN luminosity function has experi- 
enced pure luminosity evolution. It should be noted that 
we came to this tentative conclusion using a fiducial (al- 
though realistic) intrinsic AGN spectrum. In future work 
we plan to perform a more self-consistent analysis by us- 
ing INTEGRAL data on the hard X-ray spectra of nearby 
AGN. We also note that it is possible that the hard X-ray 
spectra of nearby Seyferts are in reality somewhat differ- 
ent from those of distant powerful quasars, since to a first 
approximation the intrinsic (unabsorbed) spectrum is ex- 
pected to depend on the b lack hole mass, accretion rate 
i Shakura & Sunvaevlll976l) . and spin. Observations with 
future hard X-ray telescopes will permit direct tests of 
whether the hard X-ray spectra of quasars are similar to 
those of local Seyferts or not. 
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